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ABSTRACT 


This publication presents a novel multireceiver configuration for the purpose of 
carrier arraying and/or signal arraying. Such a problem arises for example, in the 
NASA Deep Space Network where the same data-modulated signal from a spacer 
is received by a number of geographically separated antennas and the data detection 
must be efficiently perfoimed on the basis of the various received signals. The proposed 
configuration is arrived at by formulating the carrier and/or signal arraying problem 
as an optimal estimation problem. Two specific solutions are proposed. 

The first solution is to simultaneously and optimally estimate the various phase 
processes received at different receivers with coupled phase-locked loops (PLLs) wherein 
the individual PLLs acquire and track their respective receivers' phase processes, but 
axe aided by each other in an optimal manner. The coupled PLL estimator is followed 
by simple addition of the quadrature-phase components of the mixers from various 
receivers for the purpose of data detection. Such signal combining, of course, requires 
the transmission of high information rate (the rate related to the data transmission rate) 
baseband signals from all the receivers to one of the selected receivers, as expected. For 
relatively strong correlation among various phase processes, the detection performance 
as measured in terms of radio loss is near optimal; in this case the phase error processes 


themselves are also highly correlated. 

However, when the phase processes are relatively weakly correlated, coupling 
among loops is of only marginal significance (in the limit the improvement due to 
coupling is zero for zero cross-correlation). In this case, most of the gain is expected to 
result from the post-loop combining. However, simple combining of baseband signals 

iii 


is far from optimal in this case as the residual phase errors from different loops are 
highly uncorrelated. For this case, and for the case of relatively high values of symbol 
energy-to-noise spectral density ratio (F»/JV„), we propose a novel configuration for 
combining the data-modulated, loop-output signals. 

The scheme can be extended to the case of low (E s /N 0 ) case by performing the 
combining/detection process over a multisymbol period. Such a configuration results 
m the minimization of the effective radio loss at the combiner output, and thus a 
maximization of energy per bit to noise-power spectral density ratio is achieved. 
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1. INTRODUCTION 


Recently, there has been some interest in a carrier arraying system to combine 
the received carrier signals at geographically separated antennas with different carrier 
phases to improve the overall carrier tracking performance [1-3]. The expected improve- 
ment occurs both in terms of the rms phase jitter and radio loss or the carrier- to- noise 
power spectral density ratio (CNR) margin. 

The scheme of [1-3] treats one of the receivers as “master” receiver that tracks its 
received phase noise process by a comparatively large bandwidth loop. It is aided by 
other loops in combining the signals received from other receivers at an intermediate 
frequency (IF). In turn, the “master” receiver broadcasts its voltage-controlled oscilla- 
tor (VCO) output signal to other receivers so as to aid their tracking processes. Thus, 
the coupling among various receivers in this configuration occurs at both RF and IF 
stages. If the phase processes among various receivers are strongly correlated, so that 
the loop noise bandwidth of all receivers other than the “master can be reduced enough 
to achieve negligible tracking errors, then the IF combining occurs almost coherently 
and the system (first receiver) is expected to perform optimally. 

The above scheme, however, seems to have been arrived at intuitively and it is 
not clear if the performance is optimal. In fact, for the case of equal antenna size 
and receiver noise temperatures, one would expect some symmetry in the receivers 
coupling that is absent in this scheme. Moreover, as will be shown in this publication, 
if the phase processes are not very strongly correlated, then signal combining in this 
manner can result in a very significant loss of optimality, in terms of both CNR margin 
and data detection performance. In an alternative scheme presented in [4], individual 
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carrier signals are tracked by independent PLLs and the (linear) signal combining takes 
place at the baseband (subcarrier) level or at the symbol stream level just prior to the 
data detection. 

In this publication, we formulate the problem of parameter estimation for an an- 
tenna array as an optimal estimation problem in a way similar to the case of a single 
antenna [5-9]. The proposed configurations arrived at in this manner consist of two 
distinct stages as shown in Figures la and lb, wherein the first stage provides mutual 
coupling among the loops for joint estimation of various phase processes, while the 
second stage performs signal alignment and combining. While the two configurations 
are identical in terms of the first stage, they differ markedly with respect to the second 
stage. Whereas the second stage of the configuration I of Figure la involves linear 
combining of various signals, the configuration of Figure lb achieves both an explicit 
phase alignment and the signal combining. The optimization technique also yields the 
actual coefficients of coupling for these configurations. 

In practice, one could of course use the proposed structure, but the actual co- 
efficients may then be related to more intuitive parameters such as the loop noise 
bandwidths of various loops. Another advantage of the proposed scheme is that under 
relatively strong correlation among various phase processes, the mutual coupling among 
receivers is at baseband only as against at IF and RF levels in the previous scheme. In 
the proposed scheme, the coupling among receivers is via baseband error signals, and 
in the absence of coupling (simply by disabling the interconnections among various 
receivers), the system reduces to a decoupled configuration, i.e., individual receivers 
track their respective phase processes. This may be a desirable feature and is not 
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common to the scheme of [1-3]. Moreover, a consequence of this approach is that the 
system explicitly takes into account any correlation among various receivei parameters 
as noise temperatures. The scheme of [1-3] does not suggest any mechanism for this. A 
decision-directed version of the proposed scheme can be used for signals with relatively 
high E s /N 0 . It consists of the simple addition of the quadrature-phase components 
of the mixers from various receivers for the purpose of data detection. Such signal 
combining, as expected, requires the transmission of high information rate (related to 
data rate) baseband signals from all the receivers to one of the selected receivers. 

It is also to be noted that for the case of weak correlations among the phase 
processes, the scheme of [1-3] would not be satisfactory in that the phases of the IF 
signals at the input to the combiner are uncorrelated. Adding such signals would 
result in a drastic performance loss compared to the case of strong coupling. For 
this case, the proposed configuration II effectively averages out the phase errors at 
various receivers, thus performing an explicit phase alignment in addition to the implicit 
alignment achieved via mutual coupling and resulting in the same performance as for 
the strong coupling case. The scheme is applicable to the cases of both the residual 
carrier and the decision-directed (suppressed carrier) loops. Thus, the scheme provides 
a novel method for optimum signal combining under considerable phase jitters among 
estimated phase processes. 

In this technique, the high information rate (related to data rate) signals are 
transmitted at some IF frequency for the purpose of optimum nonlinear combining 
at one of the selected receivers. Alternatively, these signals may be transmitted at 
baseband and modulated at the central receiver by an IF frequency for the purpose of 
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phase aligment. Figures la and lb may be of help in the conceptual understanding 
of the two proposed schemes. In the next section, we derive these structures by the 
application of optimal estimation theory; we also present their more detailed schematic. 

Detailed performance evaluations of the proposed schemes are presented in the 
following sections of this publication. These results show that for the case of very high 
correlation among various phase processes and equal CNR among the N receivers, there 
is an improvement of lOlog^Q N dB in terms of effective CNR due to arraying, as is 
intuitively expected. However, the more important and perhaps surprising aspect of 
the presented results is that with the proposed arraying scheme, the same improvement 
can be obtained for the complete range of correlation coefficients. 

The decision-directed versions of the proposed coupled-receiver configurations in- 
corporate a transition detection scheme in the loop so as to provide a near-optimum 
estimation/detection performance even under the conditions of relatively high phase 
errors that may be encountered during acquisition. 
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2. SIGNAL MODEL AND RECEIVER STRUCTURE 

We consider the problem of estimating the phase processes 9i(k),i = 1 , 2 ,... A 
from the sampled version of the received carrier signals y'(k ), i.e., 

y\k) = AV2Sln(u c t k + 9\k)) + t7 , '(fc); i = 1,2,..., N (1) 

where is the fcth sampling time, ui c is the known carrier frequency and v W, 
observation noise, is the sampled version of a narrow-band, zero-mean, white- Gaussian 
noise process v'(t). In (1), without loss of any generality, we have assumed that the 
carrier amplitude is equal for all the receivers, while the additive noise variance may 
differ for various receivers. This is in view of the fact that the performance of the 
estimation algorithm depends only on the signal-to-noise ratio and not individually on 
the signal and noise powers. However, with minor variations, the following estimation 
algorithm includes the case where A is replaced by A*; i = 1,2, ... N in signal model 
(1). The phase processes 9'(k) at various receivers are assumed, in general, to be 
partially correlated processes (in the limit these may be completely uncorrelated or in 
the other extreme may be identical) and are modeled as 

9\k) =Hx i (k); V = [1 0...0]; i = l,2,...,JV (2) 

In (2) above, x’(fc) denotes the state vector of dimension n s associated with the ith 
phase process. For the case of a single receiver configuration, the state vector x‘(k) is 
modeled as [8]: 

x'(k + l) = $ a x i (k) + w i (k) (3) 

where is an n s x n 9 matrix and ud(fc) is a zero-mean, white-Gaussian noise process 
independent of v l (k). 
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For the multireceiver configuration considered here, we model the state dynamics 


as follows: 


‘ x x {k + 1) ' 




" w^(k) ' 

x 2 (k + 1) 

= ($ s ®F) 

x 2 (k) 

+ 

w 2 (k) 

,x N (k + 1). 


« 

sT" 

1 


I 

S 
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where <8> denotes the Kronecker matrix product (the definition is included in Appendix 
A), F is an ( N x N ) matrix, and the dynamic noise processes, w'(k), are in general 
mutually correlated. Thus, the processes x'(k) may be coupled both via the matrix F 
and via the cross covariances of w l (k), i = 1,2, . . . , N. In a simpler model, F may 
be taken to be an identity matrix, and an appropriate structure may be imposed on 
the covariance Q of w(k). The latter denotes the Nn s dimensional noise vector in (4). 
Denoting by x(k ) the composite state vector for the multireceiver configuration, one 
may write (4) in the following compact form: 

x(k + 1) = $ x(k) + w(k) (5) 

with $ = ® F. We assume that w(k) is a zero-mean, white- Gaussian process 

independent of v(k) = [{u 1 ^)} 7, . . . {v N (k)} T ] T , where T denotes the matrix transpose. 
Thus, 

F[F(fc)] = 0, E[w(k)] = 0 

( 6 ) 

E{v(k)v T (k)\ = R, E[w(k)w T (k)] = Q; ^(t)» r (i)] =0 ; j, * = 1,2,... 

The measurement equation (1) may also be written in the following vector-matrix form: 

y(k) = h(x(k), k ) + v(k ) (7) 
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where 


y(k) = [; y\k ) . . . y N (k)] T , v(k) = [^(k) . . . v N (fc)] T 


h(x(k),k) = [h 1 (x,k)...h N (x,k)] T 

( 8 ) 

h\x, k ) = Ay/ 2 Sin(w c i fc + 0‘(*O) 

= A\/2. Sin(uj c t k + F x\k^ ; i = l, 2 ,...,iV 
As shown in Appendix A, the extended Kalman filter equations for the signal 

model (5-8) are given by (Al), (A15-A17), (A5), and (A14).Under the assumption of 

a small estimation error, x(k + 1 /k) = x(k + 1 ) - x(k + 1 /fc), the filter equations may 

be simplified as follows: Representing the bandpass additive noise v'(k) in terms of its 

baseband quadrature components v\(k) and v‘ q (k) as 

v'(k) =v\(k) Sin (u c tk + £'x , (k)) + v q (k) Cos (uj c tk + £ x ( k )) 


the prediction error r) l {k + 1) given by (A17) may be approximated by 

r) l {k + 1 ) = A£ T x'(k + 1/k) + — 7 = v](k + 1 ); * = 1,2, ...,N (9) 

We note that 77 (fc + 1) = [r ? 1 (k + 1) . . . rj N {k + l)] r of (9) is also the one-step- ahead 
prediction error for the following model: 

y(k + l) = AHx(k + l) + ^=v,(k + l) (10) 


where 


■£ t 

0 ... 

... 0 ■ 

0 

F ... 

... 0 

. 0 


... F. 


F = [1 0 . . . 0 ] 
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and v,(k + l) is the vector consisting of v\ as its jth elements. It is also clear that under 
the assumption of small estimation error, the filter equations of Appendix A reduce to 


the following standard linear Kalman filter equations: 

x(k + 1 jk + 1) = $ x(k/k ) + K(k + l)r](k + 1) (11°) 

V'(k + 1) = V2y\k + 1) Cos (©*'(* + 1)) (116) 

&(k + 1) = u c i k+1 +e T x t (k + l/k) (11c) 

I<(k + 1 ) = AP(k + 1 /k)H S~\k + 1) (Hd) 

P(k + 1 jk + 1) = P(k + 1/k) - A 2 P(k + 1 /k)H T S~ l {k + l)H P(k + 1/k) (lie) 

P(k+ l/k) = $P(k/k)<f> T + Q (11/) 

S(k + 1) = H P(k + 1 jk)H T + R (ll^) 


where R is the covariance of the measurement noise vector v^k + l)/\/2 appearing in 
(10) and is equal to R/2. Note that in the derivation of the filter equations in Appendix 
A, the measurement noise covariance matrix R is assumed to be diagonal. However, in 
view of equations (9-10), the recursions (11) are applicable to the case of nondiagonal 
R as well. 
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3. OPTIMUM ESTIMATION VIA COUPLED, PHASE-LOCKED LOOPS 

The various antennas may be geographically separated in most situations, includ- 
ing space exploration. Therefore, for the purposes of proposed implementation, we 
represent the solution (11) as a set of interacting phase-locked loops with various PLLs 
associated with their respective receivers and the interaction among them occurring 
over low data rate communication links. Such an arrangement is illustrated in Fig- 
ure 2a for the case of two receivers. The receiver configuration of the figure may be 
simplified by replacing the components of the Kalman gain matrix by its steady-state 
value, which may be precomputed. In addition, the feedback gain vectors 12 and A 
(superscripts refer to the specific block in the partition of the Kalman gain matrix) 
may equivalently be replaced by scalars to bring the organization of the figure into a 
more conventional form. 

Considering first the case when F \2 = F 21 =0, and denoting by 0 (z) the Z 
transform of the input 9 l (k) to the phase modulator in Figure 2a (in actual implemen- 
tation, the phase modulator is replaced by NCO, whose input is the difference signal 
(0i (fc) - - l))(K v Ai)“ 1 with A t being the lead time and I< v equal to the NCO 

sensitivity), then O 1 ^) is given by 

e 1 ^) = PFutfsiz I - Fn*,]- 1 K n r)\z) + i T F u $s[z I - Fn^s}' 1 K 12 r] 2 (z) (12) 
It may also be expressed in the following more convenient form: 

e 1 (z) = G n (z)rj 1 (z) + G u r 1 2 {z) 

G u (z) = e T F u $ 3 [z I - Fn*.]- 1 # 11 , G 12 (z) - PFn*.[z I ~ Fu*.]- 1 * 12 (13) 
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For the case of n s = 2 and with 


$, = 


T 

1 


the filter transfer function G n (z) may be obtained as 


_ Fu(K}'+ K}' T)z-F*, K }• 

° W (14) 

Similarly, G 12 {z ) is obtained by replacing K 11 by K 12 in (14). 

This formulation ignores the internal state of the filter related to the initial condi- 
tion ^ 1 (0/0). However, by selecting the initial states of the two filters G n (z ), G 12 (z), 
such that their sum is equal to x 1 (0/0), and with the proper implementation of the 
filter, this problem is eliminated. Figure 2b shows an equivalent implementation of the 
receiver coupling for the case of F 12 = F 2 1 = 0. The filter transfer functions G 2l (z) 
*^d G (z) have a form similar to (14) for the case of n — 2. Figure 3 shows a digital 
filter implementation wherein the filter states are related to x(0/0) as discussed above. 

For the case of F i2 , F 2 i ^ 0, it is required that the filter states of the second 
receiver also be fed back to receiver 1. In this case, the implementation of a two-receiver 
configuration is depicted in Figure 4 wherein Si and S 2 are multivariable digital filters. 
A detailed schematic diagram for filter S\ is shown in Figure 5. The implementation 
of filter S 2 is very similar. 
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4. AN EVALUATION OF THE PROPOSED CONFIGURATION 

Under the assumption that the estimation errors are relatively small, as is the case 
under high CNR, the phase detector may be approximated by linear characteristics and 
the performance of the system is given by P(k + l/k) of equation (11). As k -► oo, the 
matrix P(k + 1/it) approaches a steady-state value denoted by P, which is the solution 
of the following algebraic matrix Riccati equation. 


P = $ 


p-a 2 p h t s- 1 h p $ t + q 


(15) 


As shown in [10], the solution P is given in terms of the eigenvectors of the following 
2 n dimensional matrix Hf. 

[ $- T H T R~ l H 

H = (16) 

f Q $ + Q *- t H t R-'H 

Moreover, the eigenvalues of Hf are such that the reciprocal of an eigenvalue is also an 
eigenvalue. Assuming that the eigenvalues are distinct, then Hf may be factorized as 


H f = W D W' 1 


(17) 


where D is a diagonal matrix consisting of the eigenvalues of H f along its main diagonal, 
with the first half (n) entries having their magnitude greater than one, and W being a 
nonsingular, eigenvector matrix. With the matrix W partitioned into nxn blocks as 


below, 



L W2i 

then the steady-state solution P is given by 


W 12 

W 22 


p = w 21 wr 2 a 


(18) 


(19) 


11 


The steady-state filter error covariance matrix is then obtained by substituting P for 
P(k + 1/k) on the right hand side of (lie), as 

P F = P- A 2 P P (20) 

Substitution of P for P(k + 1/k) also yields the Kalman gain matrix K from (lid, g). 

Knowledge of A then permits the implementation of the (time-invariant) digital filters 
in Figures 1-5. 

In Appendix B, we consider coupled, dual-receiver, configuration-employing type 
I and type II loops in more detail. For type II loops, which are of more interest 
here, equations (B15)-(B22) present an analytical solution in a convenient form for 
the matrix W of (18) as a function of the various parameters, rq, 7 q, T, etc., 
related to the joint phase processes model (Bll). The filter error covariance matrix 
P for a two-receiver configuration, with each receiver employing a second-order loop, 
is then evaluated on the basis of equations (B15)-(B22), (19), and (20) for various 
values of these parameters. One of the most important quantities of interest is the 
ratio of the phase error variance of a receiver operating under coupled configuration 
to that obtained without any coupling while keeping the statistics of various processes 
the same. To keep the presentation simple (where only a relatively small number of 
parameters vary), we assume that the process noise in each receiver has equal variance 
0n = 022 in (Al). This is an appropriate assumption for the present application. 
Denoting by P f* the steady-state filter error covariance matrix for the kth receiver 
operating independently, this performance measure for the kth. receiver is simply equal 
to 10 log 10 {P F i(l, 1 )/Pf( 1> 1)} and 10 log 10 {P f2(1,1)/Pf(3,3)} respectively for k ~ 


1,2, and is evaluated for various values of the parameters T, p — 512/911 0 a / a v 

( qn/T 2 r k ) with r k equal to the observation noise variance for receiver k. The notations 
a\ and a 2 v have been used in [7,8] and have been introduced here for the purpose of 
referring to the results obtained in [7,8] for the case of a single receiver. As shown in 
[8] for the case of a single receiver, the parameter c^/erj; is decisive in determining the 
optimal loop-noise bandwidth of the phase-locked loop. 

A more effective performance measure is the required increase in the received 
signal-to-noise power ratio for a single receiver, in order to achieve the same phase 
error as obtained by the dual receiver. The performance of a single receiver has been 
analyzed in some detail in [8] and some of those results are reproduced here in Figures 
6-8. Figure 10 plots the normalized, two-sided bandwidth 2 B p of the optimum phase- 
locked loop as a function of where B p = B L T , with B L denoting the actual 

loop-noise bandwidth in Hz. From these figures, it is observed that for a fixed o\, 
reduction in a 2 results in increased bandwidth, and thus some of the advantage of 
reduced a 2 is offset by increased Bl as far as the phase error variance is concerned. 
Thus, whereas in the absence of the process noise (a 2 = 0) and a fixed loop bandwidth, 
the phase error variance is inversely related to the input SNR, this is not so for the 
case of a 2 ^0. As shown in [8] for the case of a 2 ± 0, the optimum loop performance 
may be approximated by 


( 1.32(a 2 Ja 2 v ) 


Pfi(1,1)K^= < 


2\.22 


l 1-4 


2 \.25 


T = 0.1s 
T = 0.01s 


In general, for any value of T, the performance for an appropriate range of {a 2 JaD 
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may be approximated by 


Pfi(1, 1 )/<rl T = a(al/crl) b 

for some constants a, b. With p u p 2 representing the values of P Fl (l, 1) for a\ equal 
to a \ x i a v 2 respectively with cr 2 constant, then 

«/<) = 
or 

10 log 10 (alj a\ 2 ) = 10 log 10 (p!/p 2 ) X (1 - 6) _1 

Thus, the improvement in terms of the signal- to-noise ratio is (1 — 6)~ ^ times the 
improvement in phase error variance. For b equal to .25, this additional factor is about 
1.33. For a constant signal amplitude, A = 1, corresponding to the received signal 
power P 0 = 1, A 2 /a 2 = 2 P 0 T/N 0 , and thus the improvement in al is also equal to the 
improvement in terms of the carrier power to noise power spectral density ratio (CNR). 

The curves labeled “precombining” in Figures 9 and 10 depict the performance 
improvement as defined above for the case of T = 0.1s, r x = r 2 = 1.0, and for p 
= -995 and .999 respectively. Note that as the system is assumed to be linear, the 
phase variances can be normalized by the additive noise variances; the performance 
improvement depends only upon the ratio T\jv 2 , and not on the individual values of 
r 1( r 2 . Figure 11 plots the corresponding result for the case of T = 0.5s. As may 
be inferred from these figures, an improvement of about 1.7 dB in terms of phase 
error variance <r^ is obtained for the case of p — 0.999. Similar results are obtained 
for different sampling rates. Note that this information is adequate for ascertaining 
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the actual phase error variance for the coupled loops, with the performance for the 
single loop available in [7,8] and reproduced here for reference. Figure 12 plots the 
performance improvement as a function of <7 2 Ja 2 v for a fixed value of p — -.999 and T — 
0.5s. As may be observed from Figure 12, when a 2 Ja 2 v = 0.1, the improvement due to 
arraying in terms of a 2 is 1.7 dB, while the effective improvement in terms of CNR is 
2.1 dB. Performance when T = 0.1s or .01s is not significantly different. All the results 
presented in this section refer to the “precombining” results in Figures 6-32 and pertain 
to the proposed scheme 1 (ignoring the effect of a linear signal combiner) in contrast 
to the “post-combining” results also depicted in these figures. The latter performance 
curves, which pertain to the proposed scheme 2 (involving an explicit phase alignment) 
are discussed in a subsequent section. 

With the performance improvement (due to arraying) curves in Figures 9-12, along 
with the actual performance curves for a single receiver presented in Figures 6-8, one 
may determine the actual phase error variance for the arrayed system. For exam- 
ple, with (a 2 Ja 2 v ) equal to 1.0 and T = 0.1s, the phase error variance is equal to 
1.28(AT 0 /2P 0 ) without arraying, while with a dual receiver it is equal to .79(iV 0 /2P 0 )- 
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5. PERFORMANCE COMPUTATION VIA RECURSIVE SOLUTION 

The closed-form solution (16-19) of the algebraic Riccati equation, while of general 
nature, is cumbersome to compute for higher system dimensions. For a system dimen- 
sion higher than 4 (matrix H, dimension greater than 8), it is no longer feasible to 
come up with simple scalar equations as described in Appendix B for the case of a dual- 
receiver configuration, and thus one has to compute the eigenvalues and eigenvectors of 
H, by some numerical algorithms. As such computations are numerically sensitive, we 
instead evaluate the performance by a recursive solution approach. In this approach, 
one simply solves the recursions (lie, f) for k = 0, 1,2, ... M with M selected such that 
the Euclidean norm ||AP|| of the matrix P(k + M/k + M - 1)- P(k + M - 1/k+M -2), 
normalized by the norm of P(k + M- 1/k + M -2), is smaller than some specified value 
e. In these recursions, e is selected equal to 1(T 5 . The initial value P( 0/ - 1) is taken 
to be a diagonal matrix with its diagonal elements equal to 10 s . The recursions (lie, 
f) are performed for the first 50 iterations without the test, after which the recursion is 
terminated subject to the condition ||AP||/||P|| < e. Figures 13-15 plot the results in 
terms of performance improvement as a function of and for three distinct values 

of T. Comparison of Figures 12 and 13 shows that both methods yield the same result, 
thus cross-verifying the computations. 

Figures 16-19 plot the performance of a three-receiver configuration for the case 
of equal CNR for each receiver. Several values of the sampling period and <r 2 Ja 2 v have 
been considered. A few representative cases corresponding to T = 0.1, .5, .01 and 
(al/al) equal to .1 and 10 are presented here (these correspond to a loop bandwidth m 
the range of 0.2 Hz to a few Hz). For all of these cases, the performance improvement 
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over a single receiver is about 2.6 dB in terms of phase error variance and approximately 
3.4 dB in terms of CNR. 

It is also of interest to consider the case of unequal CNRs at various receivers. As 
a representative case, we consider the situation in which some of the receivers have 6 
dB lower CNR compared to the others. Such is the case within a good approximation, 
for example, with the NASA Deep Space Network with 70/70/34 or 70/34/34-meter 
antenna receiver stations. The performance for this case, with T = .01 s and = 

1.0, is presented in Figures 20-24. As is apparent from Figure 20, which depicts the first 
case, there is an improvement (for p = 0.999) of about 2.6 dB for receiver 1 in terms of 
effective CNR. When compared to the improvement obtained by a two-receiver system 
with equal CNRs, this represents an additional improvement of 0.5 dB attributed to 
the third, smaller antenna station. We also note here that as for the case of a two- 
receiver configuration, the performance improvement is not sensitive to or T. 

For comparison, we plot the performance improvement for the above case, but with 
a al°l = 0.1 and T = 0.1 in Figure 25. Comparisons of the results in Figures 20 and 
25 show insignificant differences in the two cases. For the case of unequal CNRs, the 
parameter cr a /a v is measured with reference to the receiver having the highest CNR. 

The performance improvement in the receivers with lower CNRs is even more spec- 
tacular when compared to the receiver with the highest CNR. As shown in Figure 22, for 
the above case (70/34/34-meter antenna stations), the improvement for the 34-meter 
antenna receiver is about 5 dB in terms of effective CNR. For the case of 1:1:4 ratio 
among noise variances (70/70/34-meter configuration), the performance improvement 
for the smaller station is 5.8 dB, as shown in Figure 23. Figure 24 presents the result 
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for the 1:4:16 ratio, in which case the improvement for the smallest station is 7.1 dB. 
Note that all these improvements are at p = .999 and may be a little higher or lower 
depending upon the actual value of p. 

Figures 26-32 plot the performance for the case of various four-receiver configu- 
rations. For the convenience of presentation, we fix the values of parameters T and 
al/al at .01 s and 1.0 respectively. Similar performance curves are obtained for differ- 
ent values of these parameters. As depicted in Figure 26, for the case of equal receiver 
noise variances for all the receivers, and for p = .999, one obtains an improvement of 
about 3.1 dB in terms of phase error variance and 4.1 dB in terms of effective CNR. 
Figures 27 and 28 plot the performance improvement for receiver 1 for two different 
cases of unequal CNRs among various receivers. The effective improvement for the case 
of a 1:1: 1:4 ratio among various receiver noise variances is 3.6 dB, while for the case 
of a 1: 1:4:4 ratio, the improvement is about 2.9 dB. As is also true for the case of a 
three-receiver configuration, the CNR improvement is spectacular for smaller antenna 
receivers. As is apparent from Figure 29, in a 1:1:4:4 ratio situation, the receiver with 
the smaller antenna achieves an improvement of about 6.1 dB. The corresponding im- 
provement for the case of a 1:4:4:4 configuration is about 5.4 dB. Figures 31 and 32 
plot the results on the performance improvements for the receiver with the smallest 
CNR for two more configurations with 1:1:1:4 and 1:4:4:16 ratios respectively. These 
improvements are about 6.5 dB and 7.5 dB respectively. 

It may be noted that while for the case when the same information (data) is 
carried by all of the links, the performance of only the largest antenna receiver (or 
actually the aggregate performance of all the receivers) would be of importance, there 
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are other interesting situations where the performance of all individual receivers is 
significant. This is the case, for example, with the frequency reuse (and/or polarization 
reuse) telecommunication satellites with spot beams, wherein different information is 
transmitted to different receivers but over the same carrier frequency. In this case, a 
smaller earth terminal can dramatically improve its performance (in terms of carrier 
recovery ) by proper coordination with a larger terminal. The coordination can be done 
via inexpensive, low-capacity, terrestrial links. These links need to carry only very low 
rate information representing the loop correction signals. 
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6. A NOVEL SIGNAL-COMBINING CONFIGURATION FOR 
COUPLED RECEIVERS 


For the case of most interest in this publication, wherein the carrier signals re- 
ceived at various receivers are data modulated by the same information, the system 
performance can be further improved by post-loop combining of the signals. For the 
case of a data-modulated signal, a decision-directed approach is used and the data de- 
tector becomes an integral part of the loop. Figure 33 depicts the schematic of such a 
loop for the case of a single receiver for reference. Note that the superscripts on various 
signals have been removed in that figure, and the phase detector outputs (baseband) 
are given by 


r)(k) = A Sin (S(k) + 7r D(k )) + Vi(k) 

v2 

£(fc) = A Cos(8(k ) + 7T D(k )) + 4= v q (k) 

V2 


( 21 ) 


where $ = 8(k) — 6(k ) represents the phase estimation error and {D(k)} represents the 
sampled version of the data signal, i.e., 


D{k ) = bj ; (. j - 1 )M <k< jM - 1 


( 22 ) 


In equation (22), bj represents the j( th data bit and M is the number of samples within 
one data-bit period. For the purposes of carrier phase estimation, the data detection 
could be based upon a differential approach, even for the case of coherent modulation. 
This approach has the advantage of better acquisition and tracking than is feasible 
with coherent detection. This is in view of the fact that if there is a considerable 
frequency drift in the loop (as in acquisition mode), then the coherent detector in the 
loop will provide a very high probability of error. Due to this fact, the loop may not 
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lock. However, the differential approach would provide estimates of data transitions 
with a low probability of error, provided that the phase drift (resulting due to a nonzero 
differential frequency component in the loop error) is smaller compared to the minimum 
phase separation among various symbols. For the case of BPSK modulation, the signals 
corresponding to “1” and “0” differ in their phase by 7r rad. For a more detailed and 
elaborate implementation of this concept, one may refer to [11], where the acquisition 
and tracking of a very high dynamic signal are considered. 

For the case of binary modulation, one may apply the following decision rule for 
the data transition detection: With 

j M — 1 j M — 1 

M H 2 £(*) 

k=(j-l)M *=(>— 1)M 

<j) = {(i.(i) - i.O - i)) 2 + K.0) ~Ui - 1)) 2 } 1 

the decision rule is 

bj = bj - 1 if s(j) < V T 

# 1 if s(j ) > V T 

In (24), Vt represents some appropriate threshold. Note that in the absence of noise, 
s(j) — 0 or 2^4 depending upon whether bj is equal to bj-i or not, irrespective of phase 
error 6 , such that the average phase error in the two consecutive bit periods is equal. 
Thus, Vt may be set to A in the presence of noise. Alternatively, one could compare 
|£ a (j) — £ a (j — 1)| against some threshold to detect a transition. 

Note that if the data are coherently modulated, then a separate post-loop detection 
may be applied to obtain bj coherently as shown in Figure 33. The output of this 


(23) 


(24) 
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detector would be meaningful only if the loop has acquired and the phase error 8 is 
small. 

It is clear that in the above detection scheme, there would be some signal loss 
due to nonzero 8(k). In a multiple- receiver configuration, such a possible loss can be 
minimized by averaging the phase errors before a decision is made. As the phase errors 
8\k) in various receivers corresponding to i = 1,2, . . . , N are only partially correlated, 
the average error would have a smaller variance than the individual errors 8 (k) . Figure 
34 shows a detailed schematic diagram for the scheme that achieves the desired objective 
for the case of two receivers. This has been referred to earlier as scheme 2 (Figure lb). 
In the first instance, we assume that the two receivers have equal CNR. Figure 34 
depicts the schematic diagram of an implementation that effectively averages out the 
phase estimation error 8{k ) and the noise associated with the two receivers for the 
purpose of optimal estimation. In the figure, uo d represents some appropriate frequency 
introduced for implementing the frequency (and phase) division. Representing the 
receiver input noise v'(k) (after being averaged over one symbol period) in terms of its 
baseband components v}(k) and v' q (k ) as 


v x = v\(k) Sin [u> c tk + 8 x (k) + 7r D(k)] 

+ v‘(fc) Cos [to c t k + 8 2 (k) + t r D(k )] 
then the output of the rf mixer 1 in Figure 34 is given by 


(25) 


*?b0) = 


A +-)=«.■« 


Sin [u d t k + 8 l (k) + 7 r D(k)] 


+ 4= fij(fc) Cos [u d t k + 6 l (k) + 7T D(k) } 

a/2 


(26) 


For the case of medium-to-high E 3 /Nq ratio, the output of the limiter/frequency divider 
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1 denoted rjo{k ) is approximated as 


rioW = A Sin I i(u , t t k + 0\k) + x 0(1) + «>(*;))) 

«i(*) = tan 

Similarly, the output of the lower frequency divider in Figure 34 is given by 


(27) 


?!,(<:) = ACosfi^f* + ~e(k) + *D(k) + ^(Jt)} 


(28) 


For the approximation (27) to be valid, the signal-to-noise ratio 2 A 2 / E[v }} 2 = 
2 E s /N 0 ( Es denotes data-symbol energy) must be much higher compared to 1. In 
general, for the case of MPSK signaling with M distinct phases, this ratio is given by 
(2Eb/N 0 )log 2 M. Thus, for M = 4 (corresponding to QPSK signal) and E b /N 0 — 0 
dB, the signal-to-noise ratio 2 A 2 j E[v }} 2 is equal to 6 dB and may be adequate. Higher 
values of M, of course, lead to even higher SNR. For those cases where 2E S /N 0 is not 
sufficient, the scheme may be suitably modified by extending the observation period to 
more than one symbol. The signal combiner combines the {tjd, £d} signals from the 
two receivers according to the following trigonometric identities: 

T]i(k) = r] + 7?!)(fc)££>(A:) 


= ASm{u d t k + 9(k) + 7r D(k ) + 9 n (k)} 


( 29 ) 


w) = e D mi - ri i D {kw D {k) 


= ACos{u> d t* + 9(k) -f 7T D(k ) + 9 n (k)} 
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with 


6(k) = ^(e\k) + e 2 (k)) ; e n {k) = \{e l n {k) + *’(*)) 

When Ud is an IF frequency, one may equivalently generate r]j(k) and £/(&) by 
mixers and bandpass filters. Finally, the complex-mixer in Figure 34 down-converts 

the signals rjj , £/ to baseband, resulting in its outputs given by 

r]{k) — ASin{/9(&) + 7r D(k ) + 9 n (k)} 


(30) 


t(k) = ACos{0(fc) + 7T D(k) + 9 „(fc)} 

For a high CNR case, 9 n (k) may be approximated by 

9 n (k) = tan~ 1 (v q (k)/\/2A) ; ^g(^) = ^ + ^g( fc )) 


Since {uj(fc)} and { v 2 (k )} are independent noise sequences, the variance of v q (k) is 
only one-half of the variance of either of these noise components. Thus, in addition to 
averaging out the phase errors ^(k), 6 2 (k), the scheme also effectively averages out the 
measurement noise (as an optimal detector should). Data transitions are then detected 
by substituting rj(k ), £(k) of (30) into equations (23, 24). 

The above scheme can be generalized to the case when the two receivers have 
unequal CNRs at their inputs. As shown in the following, the optimum weighting of 
the two phase error terms 9 X { k ) and 9 2 {k) would be a and (1 -a) respectively for some 
0 < a < 1, where in general a would be different than 0.5. This is achieved by using 
the frequency division by a -1 and (1 - a) -1 in the upper and lower limiter/frequency 
divider respectively. If necessary, a may be approximated by a rational number for the 
purpose of implementation. Note that this scheme affords the same weighting to the 
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measurement noise vj and vj. If a is close to r 2 /(ri +r 2 ), then this is also optimum from 
the viewpoint of minimizing the noise. In a more general case, the optimal weighting 
may be derived so as to minimize the probability of error, or equivalently, the ratio 
E\Qos 2 9{k)\/ E[vj\. For the case of equal CNRs, the optimum solution is a = 0.5. 
This “post-loop” combining can be generalized to a multiple-receiver configuration in 
a straightforward manner. 

With the phase errors 9\ and 9 2 having zero means and covariance matrix P, the 
weighted average phase error 9 = a9\ + (1 — a)9 2 has the following variance: 

E[9} 2 = a 2 P n + (1 - a) 2 P 22 + 2a(l - a)P 12 (31) 

In (31), Pi i, etc., denote the elements of the matrix, P. Setting the derivative of E[9 2 ] 

with respect to a equal to zero, the optimum weight a is given by 

a — ( -^22 ~ P12 \ 

\Pn + P 22 ~ 2Pi2/ 

For the more general case of N receivers, let 

9 = aid 1 + • • • + a^9 N (32) 

Then, the minimization problem consists of minimization of 

E[9} 2 = P a v ; a v = [ax a 2 ...a^] T (33) 

subject to the constraint 

oiyU = 1 ; u =[ll...l] r (34) 

Using the Lagrangian method of minimization, we obtain the following optimum solu- 
tion for a v : 

a v = (ii T P -1 u)“ 1 P~ 1 u (35) 
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Recall that the constraint (33) is necessary so as to ensure that the data-phase modu- 
lation in the combined signal remains equal to nD(k). 

In the more general case, where it is required to maximize the signal-to-noise ratio 
IT[Cos0] 2 / E[v 2 \, the optimization problem is more complex. However, the following 
approximate procedure is proposed for the solution. ^Ve first minimize the weighted 
sum E[8} 2 + 7 E[v 2 \ with respect to a v for any specified 7 : i.e., minimize 

a^P a v + J R a v (36) 

subject to the constraint (33). The solution as a function of 7 is then given by 

a, = {/(P + 7 J!)- 1 »)" 1 (P + 1 J!r 1 (37) 

With the expansion of Cos# in terms of its Taylor series, 

j£[Cos 2 #] = 1 - m 2 + ^m 4 — I 

o 

where m } denotes the jth moment of 0, E[9 p for j = 1 , 2 , In the general case, it is 

difficult to evaluate these moments precisely. For the purposes of this approximation, 
we may assume 6 to be Gaussian (following from a linearized approximation of the 
phase detector), and thus the optimization index may be approximated by (retaining 
only the first three terms) 

{l - a^P a v + (a„ P a v ) 2 } /(a^ R a v ) (38) 

and can be minimized with respect to 7 by an appropriate numerical technique. 

Figures 6-32 also plot the “post-combining” improvement in the phase error 
variance at the input to the coherent detector (see Figure 34). This is defined as 
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lOlog lo {P Ffc (l,l)/£[0] 2 }, where P Fk is the phase error variance for the A;th receiver 
operating by itself, and 6 is the averaged phase error given by (32). This represents the 
overall improvement in the phase error variance for the proposed decision-directed mul- 
tireceiver configuration 2, and for brevity is termed the “post-combining” improvement 
(in Figures 6-32). In these figures, we also plot the “post-combining” CNR improve- 
ment, defined as the required increase in the CNR for the single receiver, so as to reduce 
its phase error variance by the improvement factor {P p k (l,l)/ E[6] 2 } defined above. 
For the case of equal CNRs, the effective noise variance at the input to the coherent 
detector is 1/N times the noise variance for the case of a single receiver. Based on this 
information, for a given CNR, the total improvement in the ( E b /No ) ratio is given by 


(E b /N 0 )j = 10 log 10 N + lOlog lo {£[Cosr] 2 /£[Cos0] 2 } 


(39) 


For illustration, with E[9 k \ 2 = .5(rad) 2 , T = ,01s, ct 2 /ct 2 = 1.0, and a four -receiver 
configuration, Figure 26 shows an improvement in the phase error variance of 4. 5-5. 7 
dB, depending upon p (corresponding to an effective CNR improvement of 6-7 dB). 

With the approximation (based on Gaussian assumption on 6) that 

E[ Cos0] 2 £ 1 - <7 2 + a 4 ; <r 2 = E[9 2 } 


£|Cos«] 2 1 - + £ 


(40) 


we have 


E[6 2 ] = .13 ; E[Cos9 k } 2 = .75 ; £[Cos0] 2 £ .883 


(E b /No)i = 6.8dB for p small, 


and 


E[6 2 } = .177 ; E[Cos6] 2 £ .854 ; ( E b /N 0 )i = 6.6dB for p large. 
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We note that an alternative simpler scheme 1, with the signals f(fc) (see Figures 
33 and 34) simply averaged out, may result in a considerable signal loss compared 
to the proposed post-loop combining method of this publication. To illustrate the 
point, let us consider the case when p is relatively small, say 0.5, so that 6 from 
various receivers (even under coupled mode) are nearly uncorrelated. For example, 
with T = ,01s, = 1.0, and with a received signal-to-noise ratio yielding the phase 

error variance of 0.5 (rad) 2 , without coupling in a four- receiver configuration, there 
is an improvement of only 0.24 dB (see Figure 26) and the individual receiver phase 
error variance (under coupled mode) is thus E[d k ] 2 = .466 (rad) . The correlation 
coefficient between various phase errors 6 k is computed to be only .11, and thus even 
after “precombining,” various 0 k are almost independent. For simplicity, we assume 
that 6 k and hence Cos 6 k for k = 1, . . . 4 axe in fact independent random variables. 

With the approximation (40), we have 

E[Cos0 k ] 2 S* .751 

X> S «T = •<* 

* k=l 

Thus, there is a loss of 101og 10 .19 = -1.8 dB relative to the ideal case. Essentially, the 
spread on phase error variance partly offsets the improvement due to noise reduction. 
Ignoring the nonlinearities in the phase-locked loops, the improvement in terms of 
(Eb/No) compared to one receiver system is given by 

(. E h /N 0 )i = 6 + 10 log 10 (.66/.751) - 5.4 dB 

By the proposed method of post-dectection combining, however, E[8 2 ] = .1165 (rad) 2 
with E[Cos0] 2 = .897, thus resulting in an (Eb/No) improvement of 6.8 dB. Thus, 
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in this case, the additional gain resulting from the proposed post-detection nonlin- 
ear combining is 1.4 dB over the simpler linear combining method 1. As in arriving 
at this simple result, we have ignored nonlinearities (worse at higher variances) and 
made a Gaussian assumption on 6 in arriving at (40), the actual difference may be 
even higher than 1.4 dB. For higher values of p , the improvement in the phase error 
variance is attributed to both the “precombining” and “post-combining” processes, as 
is apparent from Figures 6-32. Taking into account the phase detector nonlinearities 
under relatively low SNR conditions, it is clear that the combined gain due to both 
predetection and post-detection would be actually much higher than 6.8 dB, especially 
at high values of p. This is so because the loop operating with E[0 2 ] = 0.5 (rad) 2 
(predicted on the basis of linear theory and without arraying) is expected to involve 
much higher degradation due to nonlinearities compared to a loop with E[6 2 ] = .177 
(rad) 2 (obtained with arraying), in view of the simulation results of [7,8] for the case 
of a single loop with phase detector nonlinearity. In fact, when the problems of cycle 
slipping associated with nonlinear behavior are also considered, it may not even be 
feasible to operate the loop at such high values of E[6 2 \. Thus, an additional gain of 
1.5-3 dB is expected over systems which involve only baseband signal combining (no 
carrier arraying) for a four- receiver configuration with a net gain in terms of ( Eb/N 0 ) 
in the range of 3-4 dB. Similar results would hold for other receiver configurations. 
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7. CONCLUSIONS 


We have presented two novel estimation/detection schemes for multireceiver com- 
munication systems for the coordinated estimation of the phase and frequency of data- 
modulated carrier signals. These should reduce the radio loss to a negligible value and 
thus effectively minimize the probability of data detection. 

Considering first the case of no data modulation (residual carrier) and formulating 
the problem in a state-space description, we have derived an estimator structure on 
the basis of optimal estimation theory. This configuration consists of coupled phase- 
locked loops with individual PLLs tracking the carrier signals associated with various 
receivers and aided by other loops on a mutual basis. Detailed results of analysis and 
computation have been presented to illustrate the improvement obtained due to such 
a coupling over that achieved for the case of independently operating PLLs. 

Various receiver configurations, including 2-4 receivers, have been considered for 
the cases of both equal and unequal CNRs at the receivers. We have considered a wide 
range of values for the ratio of process noise to additive noise variance (a dominant 
parameter in controlling the loop bandwidth) and a wide range of correlation coefficient 
p among the phase processes at the input to various receivers. The improvement has 
been measured both in terms of the reduction in the phase error variance (phase jitter) 
and also in terms of the extension of the CNR margin (reduction in the required CNR). 

It is apparent from the results presented in this publication, that for a relatively 
high value of p, loop arraying provides an effective CNR improvement of about 2- 
4.5 dB, depending upon the number of receivers in the configuration, for the receiver 
with the highest CNR at its input. This improvement is applicable over a very wide 
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range of loop noise bandwidths (process noise variances). For receivers with smaller 
input CNR (smaller antenna aperture), the improvement is more dramatic. For the 
case considered in this publication, wherein each receiver can have one of two possible 
CNRs (corresponding, for instance, to either a 34-meter or a 70-meter antenna), the 
maximum CNR improvement can be as high as 7 dB. 

For the case of relatively high p, the above scheme can be extended to the data- 
modulated case (suppressed carrier or data-aided residual carrier loop), by replacing 
the individual PLL by a decision-directed version. The data detection is achieved on the 
basis of a signal obtained by combining the appropriate baseband signals (quadrature 
components of the mixer) from all the loops. For relatively high values of p , it turns out 
from the analysis (as is also intuitively expected) that the residual phase errors from 
various loops are also highly correlated and thus simple combining of the baseband 
signals does not result in a significant loss of optimality. In the decision-directed version 
of the loop, we have removed the data modulation within the loop on the basis of 
differential detection or transition detection [9], as this provides a better acquisition 
performance than the usual matched filter detection. 

Apai t from being an optimum configuration, a major advantage of the proposed 
scheme with respect to an earlier scheme in the literature is that here the coupling 
among loops is only via low frequency baseband signals (loop error signals) as against 
coupling at the IF and RF levels required in the previous scheme. This feature is 
especially attractive when the various receivers are geographically separated by large 
distances as in the case with the Deep Space Network. However, it may be emphasized 
that the scheme is of a general nature and also has various other applications as well. 
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For example, in the case of frequency or polarization-reuse satellite configurations, the 
scheme can be applied for coordinated and thus improved phase estimation even though 
different carriers (channels) carry different data. 

For the case of relatively low values of p, both the above scheme and the scheme 
of [1-3] are of little use. As seen from the results presented in this publication, loop 
coupling by itself, in this case, provides very small improvement even though the esti- 
mation configuration is optimum. Thus, for p = 0.5 and a four-receiver configuration, 
one obtains less than a 1-dB improvement through coupling. In this case, it is essential 
to recognize at the outset that the overall objective is to minimize the probability of 
error or minimize the phase jitter in the combined baseband signal. An equivalent way 
to achieve this is to estimate individual phases with errors 8 k as before to obtain the 
baseband signals DCos9 k (Z) denotes data) and from these signals obtain a composite 
signal DCos9 where 9 is an optimal average of various error signals 8 k . The second 
scheme in this publication does indeed achieve such combining. From the analysis, it 
turns out that for low values of p, B k s are essentially uncorrelated and it is essential 
to average out 8 k in this manner. For the purposes of data detection, the proposed 
scheme also automatically averages out the additive noise. 

The proposed configuration II requires inputting the data signal at an intermediate 
frequency to a selected master receiver from all the remaining receivers. The mutual 
(two-way) coupling is the same as for configuration I, i.e., it involves only low frequency 
loop error signals. The performance of the proposed configuration II has also been 
evaluated under appropriate simplifying assumptions and has been presented in this 
publication. From these results, it is apparent that this scheme can provide a CNR 
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improvement of about 8 dB for a four-receiver, equal CNR configuration over the 
complete range of p. An interesting observation is that for the second scheme, the 
performance is better for lower values of p than for the higher range of p. This is 
directly opposite to the result for configuration I and may at first sight seem to be 
counter-intuitive. However, the result is easily explained on the basis of a fundamental 
averaging principle. 

Thus, the variance of the average of N random variables Xk, k = 1 (as- 

sumed to be identically distributed for simplicity) depends directly on the correlation 
coefficients among the variables \k and achieves a minimum for zero cross-correlations. 
Replacing %k by 9 ^ then justifies the results presented in this publication on an intu- 
itive basis also. It should further be noted that there is no known scheme for this case 
of low p in the published literature. 

One could, of course, use the more sophisticated configuration II even for relatively 
high values of p. At the cost of some complexity, this results in marginal improvement 
over the first, simpler scheme. 
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8. SUGGESTIONS FOR FUTURE RESEARCH 


As is generally the case with such estimators, the parameters of the proposed 
optimal estimators, e.g., noise bandwidths of various PLLs, are based on the knowledge 
of certain statistics of the received signal model as the ratio an d the correlation 

coefficient p. In practice, the unknown statistical parameters may be replaced by some 
appropriate bounds or estimates of these. This may result in some possible loss of 
optimality, which is dependent on the accuracy of these estimates. The estimates of 
( 7 ^ p may be obtained either on an offline or on a real-time basis, the latter 

procedure resulting in an adaptive optimal system. Development of such adaptive 
procedures is an important part of the suggested future research in this area. 

The performance results presented in this publication were derived on the basis 
of linear estimating theory and are applicable for a relatively high CNR case. Under 
relatively high CNR and tracking mode, the phase estimation errors are relatively small 
and the phase detectors in various PLLs may be approximated by linear characteristics, 
resulting in a linear model for the complete coupled system. However, for relatively 
low CNRs, a linear approximation may not be valid and the results presented here may 
represent only an approximation of the actual performance. For relatively low CNR, 
and/or in the acquisition mode, it is necessary to obtain results on the basis of nonlinear 
theory or by computer simulations. This forms another part of the proposed future 
research. Also, as indicated in Section 6, the proposed nonlinear phase alignment 
and combining technique is applicable when 2E s /No is much higher compared to 1. 
Investigation of a scheme that overcomes this limitation by extending the combination 
to a multisymbol period is also part of such a proposal. 
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APPENDIX A: Extended Kalman Filter Equations for the Signal 
Model (5-8) 

In the following, the extended Kalman filter equations are obtained by expanding 
the nonlinear vector measurement function h(x,k) into a Taylor series around the 
predicted estimate of and retaining only the linear and quadratic terms. These 

equations are then simplified by ignoring 2 to c t and the higher harmonic, terms piesent 
in the recursions as in [8, 12]. In fact, the development here is not much different than 
in [8]. With x(k + 1 jk) denoting the one-step ahead prediction of x(k + 1), then the 
extended Kalman filter equations for the model (5-8) are given by the following filter 
state equations: 


x(k + 1/k + 1) = $x(k/k) + M{k + \)v(k + 1) (^1) 

v(k + 1) - y{k + 1) - h(x(k + 1/k), k + 1) (A2) 

and filter gain: 

M(k + l) = P(k + l/k)h^a-\k + l) (A3) 

In (A3), P(k + 1/) denotes E[{x(k + 1) - x(k + l/k)}{x(k + 1) - x{k + l/k)} T ) 
and has the following recursion: 

P(k + 1 /k) = P(k + l/k ) - P(k + l/k)h.T(r~\k + 1 )h x P(k + 1/k) (A4) 

P(k + l/k) = $P(k/k)<S> T + Q (A5) 

c(k) = h ,P ( k/k- pwh- 1).[0 ; p(k/k- d] t Me) 
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where h x denotes the ( N x n ) gradient matrix and 

i,j = 1 ' 2 

is an N vector. 

To simplify the recursions (A1)-(A7), the gradients h x and ti xx are evaluated as 

follows: The N x n gradient matrix h x is given by 

’ Cos© 1 (A:) 0 0 

0 Cos0 2 (fc) 0 

h x = AV 2 (A8) 

0 0 0 Cos Q N (k). 

where ® denotes the Kronecker matrix product and t T = [1 0 ... 0] is a row vector of 

dimension n s . Similarly : P{kjk — 1) is the following diagonal matrix: 

h 

^2 : PWh - 1) = 

r{Sin0 1 (fc)Pj 1 } 0 •••0 0 1 

0 {Sin0 2 (&)P? 2 } 0 0 (A9) 

- As/2 

o 0 0 ••• {Sin e N (k)P£ N }. 

where P^ denotes the zth diagonal term of the matrix P ^ , and Pq is the diagonal version 

of the phase error covariance matrix, i.e., it is a submatrix of P(k/k - 1), consisting 

of elements with 1, ( n s + 1), (2 n s + 1), . . . (( N — l)n 3 + 1) as their row and column 

indices. Substitution of (A8) and (A9) in (A6) shows that 

a(k) = B(k) + R(k ) (A10) 

where B(k) is a diagonal matrix with its zth diagonal element given by 

D l \k) = 2P“ A 2 Cos 2 Q/k) + (P* < ) 2 A 2 Sin 2 0 , '(jt) 
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Denoting by P^s the diagonal matrix with its ?th diagonal element equal to ( P£)' 2 

a{k) = A 2 [P+ + + R(k)] + A 2 [P+ - \p<t>s\C M 

where Cm is a diagonal matrix with ith element equal to Cos20'(fc). 

Thus 

a~\k) = A~ 2 [I + D(k)C M ] _1 [P* + \p<t>s + R(k)] _1 (All) 

D(k) = [P+ + ip^s + R(k )] _1 [P* - \P4>s } ; R(k) = R(k)/A 2 (A12) 

Assuming R(k) to be a diagonal matrix, then all the matrices above are diagonal and 
each of the diagonal terms of the matrix [/ + D(k)C m] may be individually expanded 

into a Fourier series. Thus 

a~ 1 (k) ^ A~ 2 [P < / > + -P^s A R(k)] {Ao + A\Cm + ' ' '} (A13) 


The ith elements of the diagonal matrices Ao and A\ are given by 

1 


A, 


0 y /\ ~ [2>»( fc )] 2 


A 11 — 
•^li — 


1 = 


L —\ 

\D"(k)\ 2 j 


D"(k) \ ^1 - [D ‘ 

with Ci a diagonal matrix and its ith diagonal element given by 


C l 


D"(k) 


1 - £>”(*) 

[l + D»(fc)J 


then as in [6, 10], ignoring the high frequency terms, the recursion for P(k + \/k + 1) 


is given by 

P(Jb + l/fc + 1) = P{k + 1 /k) - P(k + 1 jk)H T A 2 


[A 2 P<f, + \a 2 P^ s + R(k + 1)] 1 CHP(k + 1 Ik) 


(A14) 
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with 

0 0 - 

h = o e o • • • o ; e = [i o ... o] 

.o e. 

Moreover, the filter error correction term is given by 

M{k + l)u(k + 1) = K(k + l)r}(k + 1) (.415) 

K(k + 1) = AP(k + 1 /fc)ff r C[,4 2 P* + \a 2 P+s + R(k + l)]" 1 (.416) 

and the zth component of the vector r)(k + 1) is given by 

•q'(k + 1) = y/2 y\k + l)Cos(0'(A: + 1)); i = l,2,...,N (A17) 

The Kronecker matrix product of any two matrices A and B is by definition the 
following (mp x nr) matrix, given in a block partitioned form: 

~ a n B .... a ln B 
a 2 \B .... a 2n B 

A® B = 

- &ml B .... CLmriB 

where A is an (m x n) matrix with denoting its ij th element for i — 1 
j = 1, 2, . . . , n, and B is any (p x r) matrix. 
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APPENDIX B: Steady State Performance Equations for a Dual-Receiver 
Communications System 


In the following, we present a closed-form expression for the performance of a 
two-receiver configuration for the two cases of first-order and second-order loop filters. 
In the first-order loop filter case, $ = H = I where I is a 2 X 2 identity matrix, 

and 

T ri 0 1 Tan Q 12 1 

(£ 1) 



ri 

0 * 


9 n 

9l2 

R = 



; Q = 




_ 0 

r 2 . 


_ 912 

922 _ 


The matrix Hf of (16) is given by 

r 1 0 


H f = 


1/ri 


0 

01 0 l/r 2 

9ll 9l2 (l + 9ll/ r l) <Zl2/ r 2 
1-912 922 912/n (1 + 922/^2) -I 


(B 2) 


With a few manipulations, the characteristic polynomial of the matrix Hf is given by 

\XI - H f \ = A 4 - (4 + t) A 3 + (6 + 2t + A)A 2 - (4 + t)X + 1 (£3) 

where | | denotes the determinant of the matrix, and 

(B 4) 


. 911922 - 9?2 _ IQI ■ t_ qn I — 

^ rir 2 nr 2 ’ r 2 


From the symmetry of (B3), it is apparent that if is an eigenvalue of Hf, then so is 
A^" 1 . In fact, |A/ — Hf \ may be factorized as 

|A I - H f \ = (A 2 - ciX + 1)(A 2 - c 2 A + 1) (£5) 

for some cj, c 2 such that each quadratic form has its roots with their product equal 
to 1. Comparison of (B3) and (B5) shows that Ci, c 2 are solutions of the following 
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equation: 


c 2 — (4 + t)c + (4 + 2 1 A) = 0 
or 

Cl, 2 = ^[(4 + t)± x/(< 2 - 4A] (£6) 

To compute the eigenvectors of H f, we obtain the matrix 

C = Cofact[A/ — H f] 

It is actually sufficient to evaluate the elements of the first row of C as functions of A. 
With some manipulations, these are given by 

C n = (A - l) 3 - *(A - l) 2 + (A - ^)(A - 1) + A 


C\2 = —(A - 1) 

r 2 


(B 7) 


C13 = ?n(A — l ) 2 


M 

r 2 


(A- 1 ) 


M 

r 2 


C'l4 — 9l2(A — l) 2 


With Aj and A 2 denoting the roots with magnitude greater than 1 of the quadratic 


forms in (B5), we have the following desired eigenvectors: 


1 

1 


1 

0 

1 


-W 12 - 


’ C11 " 

W21 


C12 


w 22 


C 12 

w n 


c 13 

1 

3 

to 


C13 

i 


0 

\—Xi 

-W42- 


-C14 . 


(B 8) 
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The steady-state error covariance matrix is then given by 


P = 


5 

to 

1 


'W n W12' 

W 4l W 42 _ 


_W 2 i W22. 


(B 9) 


Thus, the order of computations is the solution of (B6), followed by the solution of 
the roots of quadratic forms in (B5), and then the computations (B8) and (B9). The 
steady-state filter error covariance matrix is obtained by the substitution of P in (20) 
with A assumed to be 1, without any loss of generality. Alternatively, one may replace 
(B9) by the elements of the third row of C given below: 

c 31 = {(a-i) 2 -^(a-i)-^}M 


C 32 = bl 2 (A - 1) + 912 J /^1 r 2 


C,3 = {(A- 1 )’-^(A- 1 )^(A- 1 )} 
C34 = ( 2 C A — l) 2 + 51 2 (A — l)f / r i 

Numerical Example: 


Let 


ri = r 2 = 1; qn = <722 = <Zi 2 — 0.5 


(£ 10 ) 


Thus, from (B4), (B6) 

A = 0, t = 1, Ci = 3, C 2 = 2 

The root (of magnitude > 1) of ( A 2 - cj A + 1) is equal to 2.62 and that of ( A 2 - c 2 + 1) is 
1.0. (This is the limiting case as Q is a singular matrix.) From (B10), the corresponding 
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eigenvectors (after possible scaling) are given by 


[1 1 1.62 1.62] r ; 


[1 - 1 0 0]' 


From (B9), the steady-state covariance matrix P is computed as 


P = 


'0 

1.62' 


1 

r 

-l 

'.81 

.81' 

_0 

1.62 


-1 

i 


.81 

.8! 


The filter error covariance matrix P F computed from (20) has all its entries equal to 
.31. Without arraying, the corresponding value is computed from (15), (20) and is 
equal to 0.5. Thus, arraying results in an improvement of about 2 dB in the phase 
error variance. 

In the second-order loop filter case, the various matrices are given by 



1 

*e« 

Co 

i 

o 


1 

•s 

i 

il 

-o- 



; *.= 



O 

i 



1 

r_| 

o 
1 



'1 

0 

0 

O' 


r i 

0 ' 

H = 





; R = 




0 

0 

1 

0 


0 

r 2. 


Q = 


<hiQs Q 12 Qs 


T 2 /3 T/2 


; Qs = 


_Ql2 Qa 922^3 


_ T/2 1 


(£11) 


Substitution of these matrices into (16) yields Hf, which after simple row operations is 
modified into matrix H fm . The modifications consist of adding TxSRj to T x 3? 3 
to — y x to 3?5 and -y X to !J ?7 with denoting the i th row of Hf. The 
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corresponding matrix (A I - H f m ) is given below. 

A I-H fm = 


(A - 1) 

0 

0 

0 

-1/rx 

0 

0 

0 - 


AT 

(A — 1) 

0 

0 

0 

0 

0 

0 


0 

0 

(A — 1) 

0 

0 

0 

- l/r 2 

0 


0 

0 

XT 

(A — 1) 

0 

0 

0 

0 

(B 12) 

0 

T 

-6911 

0 

“?9i2 

(A - 1) 

A 

0 

0 


T 

79ii 

~9n 

T 

7912 

— 912 

71 £n 

2 r t 

(A - 1) 

T q\2 
2 r 2 

0 


0 

T „ 

“7912 

0 

T 

“7922 

0 

0 

(A-l) 

A 


T 

jQ 12 

~<7 12 

T 

7922 

— 22 

T g!2 

2 r ! 

0 

T <722 
2 r 2 

(A-l). 



A = ~ f(A-l)-T 

With the application of the following matrix identity [13], 

A C, 

' = \A\- \D — CA~ 1 B\ 

B D 


(B 13) 


and after a few tedious manipulations, one obtains the required characteristic polyno- 


mial in terms of x = A — 1: 

|A/ — Hf \ = a(x + l) 2 (z 2 + 6a: + 6) 2 + bx 4 (x + l)(ar 2 -f 6a: + 6) + x 8 

AT 4 9 T 2 <7] j floo 

a ~ ’ A = 9ll?22-g 12 ; 6 = — ( 1 ) 

oorir 2 6 rj r 2 

Alternatively, one may express \\I — H f \ in terms of A as 
P{ A) = A 8 + (b - 8)A 7 + (a + 28)A 6 + (8a -9b- 56)A 5 


(B 14) 


(B 15) 


+ (18a + 166 + 70)A 4 + (8a - 96 - 56)A 3 + (a + 28)A 2 + (6 - 8)A + 1 
From the symmetry of the coefficients of the polynomial P( A), it is observed that if Ax 

is the root of P( A) then so is 1/Ax for any A^ Thus 

4 

T(A) = II( A2 + C > A + !) (2? 16) 

1=1 
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where denotes product. Comparison of (B15) and (B16) shows that 

4 

d\ = ^ C{ = 6 — 8 

i=l 

E c ' c i = ° + 24 

i = 1 j=i + 1 

^ E E E C * C J Cfc = 8a ~ 12fe _ 32 

t=l j=i-fl k— j 1 

d A = C 1 C 2 C 3 C 4 = 16(a + 6 + 1 ) 


(517) 


and that ci, c 2 , c 3 and c 4 are negatives of the roots of the following equation: 


y 4 + d\ y 3 + di y 2 + d 3 y + d± — 0 


( 518 ) 


The computation of the required eigenvalues (with magnitude > 1 ) of if/ proceeds as 
follows: Evaluate coefficients df s from (A14) and (A17), and factorize the polynomial 

(A18) into the product fli=i(y + c »')- For eacb value of i “ !> 2 > 3 ’ 4 ’ obtain tbe roots of 
the quadratic equation, A 2 + CjA + l = 0, among which only the one with its magnitude 


> 1 is retained. 

With the application of Schur’s identity [13, 14], 


’A 

B' 

-1 

(A^ 1 + A~ l BHCA~ l ) 

-A~ l BH 

C 

D 


—HCA~ l 

(D -CA~ l B)-\ 


(B 19) 


and with some algebraic manipulations, the elements of the last column of the matrix 
Adj(\I — Hf) are obtained. In the computation of the Adj(XI — Hf m ) from Schur’s 
identity, we simply cancel out |A I - H fm | from the elements of (A I - H fm ) _1 - Also, to 
obtain Adj(XI-Hf) from Adj(XI - H fm ), we post-multiply the latter by the following 
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matrix: 





■ l 

0 

0 

O' 


■l 

T/3 

0 

0 - 

E, 

0 ' 

; E x = 

-T 

1 

0 

0 

; E 2 = 

0 

1 

0 

0 

0 

E 2 


0 

0 

1 

0 


0 

0 

1 

t 3 




. 0 

0 

-T 

1 . 


.0 

0 

0 

1 . 


T/3 


The elements of the last column of Adj(\I - H f ), (any other column could be selected 


as well) with Cj denoting the j th elements, have the following expressions in terms of 


= A — 1: 


C\ = 


—xB(x) 


gi 2 T 3 

187-^2 


C 2 = (x + l)B(x) 


q 12T 4 
18r!r 2 


C 3 = x(2x + 3)U(x)- — ; C 4 = — (2x 2 + 5x + 3)U(x)— — 

oV 2 O' 2 

(B 20) 

r, Jl jjf ^ n 

C 5 - -X B(x)- ; C 6 - X C(x) ^ 

C 7 = x 2 (2x + 3)t/(x)| i C 8 =x 3 C(x)-(x+l) 2 Y(x)^A 

where B,C,U,Y are the polynomials in x as follows. 

B(x) = (2x 2 + 5x + 3)(x 2 + 6 + 6) ; C(x ) = (2x + 3)(x 2 + 3x + 2) 

T 2 AT 2 

tf(x) = x 4 + (x + l)(x 2 + 6x + 6)^— ; Y(x) = x 4 + (x + l)(x 2 + 6x + 6 )^T^ 

(521) 

The eigenvectors of the matrix Hf corresponding to A;, * = 1, . . . , 4 are obtained 
by evaluating (B17), (B18) at x = Xi — 1, i.e., 


Wi = [Ci c 2 c 3 ... c 8 ] t | i=a ._ 1 


(B 22) 


The steady-state error covariance matrices P and Pp are then evaluated from (18)-(20) 


as for the case of a first-order filter. 
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Fig. la. Proposed Estimation-Detection Scheme 1. 







Fig. lb. Proposed Estimation-Detection Scheme 2. 
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Fig. 2a. Multireceiver Configuration With N = 2 (Coupling Consists of Vector Baseband Signals). 
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Fig. 2b. An Equivalent Implementation for the Case of Diagonal 

Matrix F and n = 2 (Coupling Consists of Baseband Scalar Signals). 
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Fig. 3. Digital Filter Implementation of G n (z); Filter States Correspond to x 1 (k/k). 
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Fig. 6. Normalized Loop Noise Bandwidth of the Second-Order Filter 
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Fig. 9. Improvement Due to Arraying; p = 0.995 (Close-form Solution). 
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Fig. 10. Improvement Due to Arraying; p — 0.999 (Close-form Solution). 
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Fig. 12. Improvement Due to Arraying; T — 0.5, p — 0.999 (Close-form Solution) 
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Fig, 13. Improvement Due to Arraying With a Two- Receiver Configuration; T — 0.5 (Recursive Solution) 
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Fig. 16. Improvement Due to Arraying With a Three-Receiver Configuration (Equal CNR; <J 2 a /crl = 0.1). 
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CORRELATION COEFFICIENT p 

Fig. 30. Improvement Due to Arraying With a Four-Receiver Configuration and a 1:4:4:4 Ratio (Receiver 1). 
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Fig, 34. Decision-Directed and Coupled Receiver Configuration. 
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